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ON THE POWER OUTPUT OF SOME IDEALIZED SOURCE CONFIGURATIONS WITH 


ONE OR MORE CHARACTERISTIC DIMENSIONS 

by 

Harold Levine 


§1. Introduction 

Any roster of theoretical sound wave generators includes the discrete 
point source linear array, a continuous finite-size linear aggregate, the 
plane circulat piston or annulus and a plane rectangular piston; for the far 
field directivity pattern of all these configurations, as well as their 
individual power outputs (a function of the characteristic dimension(s) and 
wave length, in time-periodic regimes), can be accurately determined if each 
source distribution possesses a uniform strength and the pistons are flush 
with a rigid baffle. The abstract nature of source descriptions which post- 
ulate a uniformity of motion (e.g., a surface or volume displacement) is evident, 
although attempts to impose more realistic boundary conditions (while retaining 
the fiction of an ideally rigid vibrating surface) encounter considerable math- 
ematical difficulties. The radiation problem becomes well nigh intractable, 
moreover, for real material surfaces whose actual motions cannot be assigned 
in advance. Such difficulties are left out of consideration in what follows, 
where the aim is merely to extend and refine the prior analysis of some model 
sources. 

The calculation of power output from a (finite) linear array of equi- 
distant point sources is investigated in §2, with allowance for a relative 
phase shift and particular focus on the circumstances of small/large individual 
source separation. A key role is played herein by the estimates found for a 
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twin-parameter definite integral that involves the Fejer kernel functions, 
viz. 


0 . 


= ] C 
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Z 
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where N denotes a (positive) integer; these results also permit a quanti- 
tative accounting of energy partition between the principal and secondary 
lobes of the array pattern. Succeeding sections §3, 4 are concerned with 
continuously distributed sources along a finite line segment or an open-ended 
circular cylindrical shell; and estimates for the relatively lower output in 
the latter configuration are made explicit when the shell radius is small 
compared to the wave length. The last section, §5, furthers a systematic 
reduction of divers integrals which characterize the energy output from 
specific line and strip sources. 
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§2. A Linear Array of Point Sources 

Consider an ensemble of N point sources which are uniformly spaced, 
at the distance d apart, along a straight line; and assume that the corre- 
sponding excitation potentials referred 
to coordinate origins at the individual 
source locations. 
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Involve a common strength (the maximum rate of volume displacement) . 

The time average radiation intensity for each isolated source is then given by 

0 . 


Vt 


(where p specifies the equilibrium density of the ambient medium) and the 
intensity for the (non-interacting) aggregate has the representation 




M 








3 = 3 , 

' •<■=<, »MJ. 


( 1 ) 


where the factors 


A| = i/ A^=e . ^ 3 =e 


, ^=« 




( 2 ) 
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incorporate a fixed phase advance oi kd between neighboring and consecutive 
sources. After the summation in (1) is effected, it follows that 



and thus the total radiated power assumes the form 



Thiessen and Embleton (1958) cite the latter measure in connection with a 

proposal to simulate multiple siren noise from suction rolls at a paper mill 

in terms of a phase-advancing linear source array. An approximate integration 

scheme suggested by them employs a subdivision of the range and an average 

2 

value for the rapidly varying function sin Nx , n >> 1 , yielding the estimate 

P/c = 3.0^M 

7T 

on the proviso that a does "not deviate by more than 1 - from the 

values 2m7r/kd , where m is an integer." In point of fact, the numerical 
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coefficient 3.07 should be replaced by it and the remaining terms of (5) 
have a questionable status. 

To study the above (power) integral without an initial manner of ap- 
proximation, it is observed that 
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through an integration by parts. On next utilizing the relation (familiar 
from the convergence theory of Fourier series) 
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an exact result follows, namely 
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The latter permits a ready determination of I if N and the number of 
terras in the sum is small; when N is large, on the other hand, the re- 
arrangement 
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may be noted, in view of the well-known Fourier series expansion 


CO 



v\-l 


whose sura is an odd (27T-periodic) function, namely 


F(^) = ; < S < 


( 10 ) 


and the remainder estimate 
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v \ = 1^+1 
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= |--Sc{n(5-Iw,x)} , (11) 
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that involves the sine integral function 


Si (x) = j 


SVM t 


<ix 


(12) 


In accordance with (7) - (11) , 
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Nci 


I(a;M)= Nci + M ^ ^ 


\A = | 
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= N Sc (iNa) - 
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; N » a « 1 , 


(13) 


There is a ready usage for the preceding result in connection with 
the directivity gain of the considered source array, (cl) , defined as the 
ratio of the maximum far field intensity (at the angle 9 = cos~^ a from 
its reference line) to the total radiated power/Air ; thus 


QU) 
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and 
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for the specific phase parameter value a = 1 
passage to the limits N , d 0 wherein 
function 




, on the basis of (13) • A 
Nd ^ £ yields the gain 



(16) 


which, as verified in §3, befits a uniform line source distribution with a 
continuous linear phase variation (kx) along its length (5,) • If a = 0 
the corresponding gain function is secured through the mere replacement of k 
by k/2 and its limiting version is appropriate to a continuous line of 
in-phase sources with overall span £ • 

Analogous estimates for values of the phase parameter a between 0 
and 1 may be given in the circumstance that kd « 1 or 

and 

+ = ^M7V + .^TV, 5«|) 


in particular, 
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= Mn - M |iT- - Sc 

where the first term predominates when the arguments of the sine integrals 
are large enough to justify the estimate Si(x)= tt/2 , x » 1 , The gain 
function inferred from (14) and (17) in the continuum limit d 0 , Nd , 
viz. 





(18) 
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extends the prior result (16) . 

Since a sharp estimate is lacking for the trigonometric sum that appears 
in the representation (7) of the (power) integral (6) an alternative analysis of 
the latter becomes warranted. A relation from the theory of Fejer sums of 
Fourier series has utility in this regard, namely 




( 19 ) 



as do the modified versions 


I M) = I (w; N) = N 

and 

) n) } = ('^+|^)r('ln;Nl) = ('i'M+i)KM 

for integer values of m • If, then, 

0 <o^< TV 


the rearrangement 
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proves advantageous, inasmuch as the first two integrals are directly re- 
ducible while the third can be suitably expanded in a straightforward manner 


for large values of N ; thus 

“• U Nk x’- “■ 
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and 
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which yields the representation 
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(22) 


Evidently the initial pair of terms in (.22) duplicates those which make up 

the prior estimate (13), while the remaining terms are collectively small 

2 1 

if N » 1 and a « 1 . Because the function esc x j analytic 

X 

on 0_<x_<a<7T it is permissible to Integrate by parts any number of 
times in the last integral, thereby generating an expansion with reciprocal 
powers of N , i.e., 




etc. 



- 12 - 


Suppose, next, that 

o, = + b ^ I bl < TT ^ 

then, recalling (20), 

I(gl)K) = 2 hK K + I(b;Ni) 


(24) 


where the second term, I(b;N)^ is describable in the manner previously 
detailed. Thus, given the configurational parameter 

- 2vin 0<<^<7v 

an estimate 

_ I Cxrv^ 

I(U;N) = J 
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M-M 




follows; and the related gain in the limit kd >> 1 (n » 1) , viz. 


(^(0- 


I(fedL;Ni) 




corresponds, as it should, to independent contributions from each of the 
N sources. 

Once the integral (6) has been suitably determined (as a function 
of the limit a) it becomes possible to specify both the whole power output 
as well as the lesser amounts radiated into different angular subdomains of 
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space, For instance, if there is no relative phase advance between sources 
of the array (0^= 0 in (4)) the output can be written as 



(where c - co/k designates the speed of sound) and thus, recalling the 


particular integrals (20) , 

P = f ■ 
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V\ = . . . 

(25) 

When 

fed 

II 
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<5<tt 



and use is made of (22) there obtains 


I(wk-|-)M) = - I(|-)N) 

= unNl - N Sc(s)^) + + 

N»| 

together with the total power estimate 

P = f — — ® hkN-NSi(M^) Lj1+... 
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Similarly, when 


“fed = TV 


it follows that 






= vi7\ Ni Ni ?c|N(n-<5)} "" — + • • N» I 


and 
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With the current stipulation regarding source phases (a = 0) the 
array is termed a broadside one, for (cf. (3)) its maximal power output 
appears at right angles to the line of support (6 = 7r/2) ; and since the 
null directions are specified by the relation 


C02.9 = ±. 





( 26 ) 


the angular range of the major lobe (m.l.) is confined within the sector 

i M?iJL ^ NiU 
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Accordingly, the energy radiated into the major lobe has the magnitude 
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which decreases as kd becomes larger. If kd - tt , i.e., the individual 
source separation equals a half wave length, the comparative ratio 


^=iSc(x,) N»l 


has a value slightly less than unity, inasmuch as the asymptotic development 
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implies that 




When kd = 2tt and the wave length is precisely equal to the separation 
between neighboring sources of the array a reduced value for the ratio 


(26), namely 
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(29) 


makes clear the existence of significant radiation outside the major lobe. 
Secondary maxima are found by means of the condition 




(30) 


,Nkd 


(that corresponds with unit absolute magnitude of the factor sin( ^ cos 9) 
in (3)), and for large N the particular determinations 6 = 0,7T (p = N) 
follows. The null directions closest to the line of the array satisfy the 


equation 




which follows from (26) on choosing p = N and yields the explicit values 
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and 


3 =T^-n/N 
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respectively. Thus, the total energy radiated (equally) into the pair of 
secondary lobes whose axial directions are aligned with that of the array 
has the magnitude 


P 5 I. = ? 




= o Sc(im) + • • • • N»l. 

^ Ur \.T^ J ' 
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and the ratio 


hi. 


p 




= Sc (in ) 


(31) 


implies (after recalling (29)) a like contribution to the output from the 
major and secondary lobes. 

An array whose radiated intensity peaks along the source line (0 = 0) 
is termed end-on or end-fire; and, since the direction of maximum far field 
amplitude is given by cos 6 = a , it appears that the specific value a = 1 
(with relative phase advance kd between neighboring sources) characterizes 
such an array. The half-angle of the major lobe can be deduced from the relation 




with the result that 
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and thus the power which leaves the array in this angular range equals 
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A simple calculation establishes 

cJ-qI I (^<>-; M ) 


as the total output from this array, whence 


Pwi. ^ N5c(rn) 
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2 1 1 

with the particular values — Si(27r) , — Si(2Tr) , and -;r- Si(27r) if kd = tt/2 

TT 7T 2tT 

TT and 2 tt , respectively. Evidently, the relative importance of the major 
lobe radiation declines after the distance between adjacent sources exceeds 
a greater wave length; by contrast, a half wave length provides the correspond 
ing reference scale for the broadside array. 
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§3. A Continuous Line Source Distribution 

If simple sources are continuously distributed along a finite line 
segment, with uniform magnitude and a linear phase variation, according 
to the density function 


= Rt,{ 


= Rc{<l,^e 

the complex space factor (})(_r) of the resultant velocity potential sat- 
isfies an inhomogeneous equation 


= 0 , 

and the companion acoustic (over) pressure p(r^) is given by 

j>(r) =-iyco4>(t;. 


The outgoing wave solution of (33), 








J ('X - nC)\ j 


_ % 








( 35 ) 



- 20 - 


has an asymptotic form 
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which underlies expressions for both the far field intensity S and the 


power output P , namely 
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and 


TK 

P = 

O 



(37) 


(38) 


Thiessen (1955) cited the latter integral in advance of the discrete 
point source study referred to earlier and, after invoking the substitution 
sin X = /ttx/2 along with an indefinite Bessel function integral 

formula, gave an exact infinite series expansion for the output which ”is 
not an attractive one to use.” However, as soon perceived by Burgess (1956), 
there is an alternative representation 
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[easily found on writing = d(- and effecting a partial integration] 

X 

that readily lends itself to evaluation or estimation for diverse magnitudes 
of the parameters k£ and a . 

A typically difficult integration over all spatial directions of the 
far field intensity need not be the sole basis for calculating the radiated 
power P ; another procedure (whose counterpart in antenna theory is better 
known) makes use of an expression for the (local) rate of energy input at 
the source itself, viz. 




(40) 


which has patent advantages if the domain of q(_r) is compact. Thus, (32), 
(35), and (40) imply that 
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and, through the introduction of a new pair of integration variables. 
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there obtains 
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with the anticipated outcome (39) after performing the final (t) integration. 

The directivity gain of a line source has previously (§2) been inferred 
from that of its discrete analogue; and the given result (18) does indeed 
emerge when the expression (39) is substituted into the formula 



Taking note of the property 

Sc (~7S ) = - St 

and of the asymptotic estimate (28) for Si(x) , the deduction from (39), 
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iniplies a lessened radiation with increasing values of the phase parameter ct • 
It needs to be observed that neither Thiessen (via display curves for P 
vs. at different values of a ) nor Burgess (in analytical work) have 

brought out the fact that a trigonometric function depending on (XkH 
enters into the power estimate appropriate to a rapidly varying phase regime. 

Side lobes in the directivity patterns are manifest for waves whose 


length is relatively small compared to the overall scale of the continuous 
line source, whereas longer waves (compared with the spacing between discrete 
source elements) can be admitted when many of the latter are present. 

If the phase variation along the source line has a centered quadratic 
behavior, as is made explicit by the inhomogeneous term 
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for (33) , the corresponding far field velocity potential 
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contains an integral which can be simplified on completing the square in 
the exponential function; thus 

ftJt cold 
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with the customary definitions 
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= -C(-9^) 
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of the Fresnel integrals. 

The net power output, obtained through an integral over all directions 
in the far field, 

p=-^^"U {inJ<N4>r'-3a3} 


is specified by 
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It is a simple matter, on utilizing the asymptotic estimates. 
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and 
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to deduce from (45) that 



so long as k£ does not assume very small values. The amount of radiation 
is thus diminished when the magnitude of the phase parameter 3 increases, 
and the occurrence of the first reciprocal power of 3 in the leading term 
of (47) may be contrasted with the second reciprocal power of a in the 
prior case of a linear phase variation along the source line. A slightly 
more intricate handling of (45) is called for in order to confirm the limiting 
value of the power when 3 0 \ once again use is made of (46) and it appears 


that 
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whence 
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in accordance with the result yielded by (39) for a = 0 . 

On adopting (40) as the means of power specification there obtains 
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as an alternative to (45) • Estimates are readily secured from (48) when 
3 0 and k£ is finite, through series expansions of trigonometric func- 

tions, and the consequences 
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(49) 


reveal a lessening of the power output associated with small phase differences 
between neighboring points of the source distribution. The analogous feature 
is also present in the case of linear phase variations specified by the para- 
meter a , where the deductions from (39) comprise 
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P(ot}“ P( 0 ) =- SwMJ + 0(o<'*'^ 


a 0 , kJl finite (50) 


and 




p (=<) - P(o) = - o< ^ 0, 

A preliminary transformation of the power formula (48) becomes 
appropriate if estimates suitable to large values of 3 are sought; and 
the first pertinent relation is given by the integral 
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Further, on designating 


(51) 


it follows that 
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permit a sequential refinement of estimates for (52) when \x » 1 ; and this 

leads, via (48) and (51), to an output characterization which isolates the 

2 

asymptotically predominant term, /83 > 6 > and contains only 

first powers of the Fresnel integrals [compare (45)]. The explicit construct- 
ion of a development for P that involves reciprocal powers of Bk£ is, 
nonetheless, a technically laborious matter. 
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§ 4. A Cylindrical Shell Source 

Consider a circular cylindrical shell with negligible wall thickness 
and open ends, whose two characteristic dimensions are A (the length) and 
R (the radius), respectively; and suppose that its curved surface is the 
site of a simple uniform source layer, specified by the density function 
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^ g«-R) 

2.7V 



0 < 9 < Zt\ 


(53) 


in terms of cylindrical coordinates r , 0 , z with origin at the geometrical 
center of the shell. The velocity potential for this source aggregate 
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possesses an evident symmetry around the axis or center line of the shell and 
leads, when combined with (53) in (40), to the power representation 
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A transformation of integration variables (cf. p. 22) yields 

tk I 

]■ = 2.^ G(|^l){i-5}A^ 

-l/t ° 


and thus 


P = ^ i ■is I (l-U I ==-^ i^ 

o o i V^+ 3/2.)^ 






-33- 


it follows that the output of the source can be recast in the form 
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The latter expression reduces in the limit R -> 0 , to 
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r f) 1 


2 . 


as befits a line source and, therefore, agrees with an appropriate special- 
ization (a = 0) of (39). 

A power estimate is readily forthcoming from (58) if 


« 1 , 
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i.e., if the radius of the shell is small compared to the acoustic wave 
length and also the length of the shell. It turns out, by straightforward 
analysis, that 
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where the first correction due to the finite radius of the shell involves 

2 

the product of (kR) and a function of k£ , viz. 


F(U) ={;SaM) + 
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with the (positive) limits 


and 
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Since the derivative of F , namely 


^ _ J 


does not vanish on the interval x > 0 the correction terms in (60) for a 
small but finite radius have a combined negative value (regardless of the 
magnitude of kl ) y and thus the shell radiates less power than a line of 
the same length* 

ilcctz 

After introducing a linear phase factor e into the source 

function (53) the power formula can be expressed in the versions 
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where the second is consequent to a reliance on the integral representation 
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Once the integration relative to the variables z , z' are effected and 
account taken of the complementary specifications 

, ?>k 
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implicit in (62), it follows that 
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and the occurrence of finite limits enables a straightforward reduction of 
the latter formula, when the inequality kR « 1 holds, through a power series 
development of the Bessel function. Thus, on employing the first two terms 
of the series, i.e.. 


= I - it 
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the difference between shell and linear source radiated powers, to the lowest 
order in the radius of the former, is expressed by 
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shell line 
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without restriction on the magnitude of . Given the particular value a = 0 , 
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it turns out that 
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in complete agreement with (60); thus, the shell with a uniform source density 
radiates a smaller amount of energy than a line with the same density, for all 
magnitudes of the ratio, length of shell/acoustic wave length. 

If the phase parameter a equals unity, then 
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where the dependence on kJl is contained in a function 
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which relates to the cosine integral, Ci x , along with the Euler constant 
Y = 0.5772 ... , and has the limiting behaviors 
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Since 





7^>0 


the inference (6P) ^ < 0 may be drawn, and thus a diminished level of radi- 

a=l 

ation is linked, once again, with finite source radius. The respective pro- 
portionalities of and (5 P)^_q to log k£ and k£ , when kS, » 1 , 

are noteworthy as is also the comparative power ratio 




p 


o (=0 




M» i 

I 
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and, in the limit R ■> 0 , 
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as is fully consistent with (41) . 
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§5. Divers Integrals 

Although the power and gain of linear and planar source models have 
long been a subject for analysis in both the acoustical and electromagnetic 
contexts, additional details and improvements on specific results are feas- 
ible, as shown in the sequel. There is a direct acoustical counterpart of 
an investigation by Clemmow (1966) on antenna gain and supergain, for example, 
wherein the resulting integrals can be displayed in a more useful manner; thus, 
a two-dimensional electric field polarized in the z-direction, say, which is 
generated by the unidirectional surface current density in the plane y = 0 , 
namely 

) =(o,o^jw» 


where 



jo, I'^l < 1 

0, hi > Mil 



constant 


( 66 ) 


corresponds with a two-dimensional acoustic field produced by a rigid strip 
radiator with the normal (y-component) velocity distribution 
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in the plane y = 0 . 
If 


constant (67) 


G(5) = n 


power radiated per unit angle in direction 6 


total power radiated 



-40- 


defines the gain for radiation into the half-space 0 < 0 < ir 
-CO < z < °° , straightforward analysis reveals that 



y > 0 


( 68 ) 


for either of the source distributions (66) or (67) , k denoting the (ap- 
propriate) wave number. In accordance with the relation 



(69) 


involving Bessel functions of zero and first orders, an expression for the 
gain (which Clemmow does not state) follows (compare (16)), 

G(y) - M ’ 

O 

that readily furnishes the respective estimates 


= I, M« 

(71) 

iiL, fci»i 

I 

Hence, a uniform current or normally vibrating strip produces a gain in the 
direction transverse to the source plane which rises in proportion to the 
strip width, provided that the latter exceeds a wave length. An alternative 
to large scale configurations for surpassing this ’normal’ gain, Gq , has 



practical merit; and this contemplates variable source distributions acting on 
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small excitation areas. Consider, in a simplified though illustrative sense, 


the feasibility of supergain G > Gq with a non-uniform strip velocity 

; \-x\<lk ] 

”N » (72) 

0 , 1-)(| > ill ^ 

in place of (67); here the v are constants and v = v . The appertain- 

n -n n 

ing velocity potential in the radiation field has a cylindrical wave form, i.e., 
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with the angular factor 
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Since 




does not depend on the parameters ‘ selection which minimizes 
the total radiated power also maximizes the gain in the direction tt/ 2 (normal 
to the strip) . Clemmow makes a choice 
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and observes that the corresponding gain 




= -I 


\ 1 — r:: 1 1 + "7: tt; 


o ^ 
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can be rendered a maximum by a suitable (real) determination of n ; and this 
yields 


'-I- 


K sw\^usib) ' ^ 
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(73) 


where is the normal gain (70) • 

Estimates for the integrals appearing in (73) are easily secured when 
k.5, « 1 and reveal that 
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to the lowest order of approximation. Clemmow also arrived at the complementary 
result A ^ 

^ I W-“) 
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(74) 
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Cwhich suggests only a small fractional increase of gain over the normal 
amount) through function theoretic reasoning. It is a more exacting task to 
obtain fuller details of the (asymptotic) expansions for all the integrals, 
when kil » 1 , by complex variable methods; hence, alternative representations 
of the integrals, which permit the requisite expansions, may be noted here and 
attention drawn to the circumstance that similar integrals are encountered in 
other problems of sound radiation. 

One of the three separate integrals in (73), connected with the normal 
gain (.68), has previously been given an alternative form, (69), which enables 
systematic estimation at both small and large values of the parameter k£ • 

The remaining pair, say 




031^5 5wi*'(^C0tS) 


enter jointly into the relation 
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(76) 


and may be rewritten as 
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and 
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(78) 


Define 


such that 


then 
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and so 

which provides the sought-after result 


( 79 ) 
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Having regard, similarly, for another two-variable function 


In ^ o 








which obeys the differential equation 
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it follows, after successive integrations consistent with the conditions 
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and 
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(80) 


Expansions of the functions , and F^ which are useful for 

small values of 3 can be secured from the respective integrals (69) , (79) , 
and (80) on substitution therein of the Bessel function power series. By 
virtue of the relations 
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and 
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it is expedient, for large values of 3 » to employ the representations 
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together with as 3 rmptotic developments of the Bessel functions, viz. 



The estimates 



and 
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follow and their leading terms jointly yield another. 


■ Go _ 


(F,(f)) 
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which confirms (74) • 

The integrals 
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(81) 


5vM even 

v \ : 

CCTL odd 




(82) 


occur in a study (Gomperts, 1974) of sound radiation from a rectangular plate 
whose vibration pattern is one-dimensional and parallel to its free edges, the 
other pair being hinged; evidently 

I -H- 


and 


F,,. ^ ^ A) = 

I 3 


F3 (^,i} = [k {F,(iv^;+nYF^Cif^)} 


while the reductions of (81), (82) for other integral values of n can be 
achieved by the means employed in connection with » ^2 ’ Gomperts calls 
attention to the * critical frequency' integral 



and to an estimate 



N) »1 


( 84 ) 


whose derivation is not given; a complete evaluation of this integral com- 
mences with the alternative version 


-4 * Jim 
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and proceeds with successive integration by parts, viz, 
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T Jp^T-Kiit; IV J,(ts)n; 


The estimate (84) is directly inferred from the first (and predominant) term 
of (85) on utilizing the Bessel function approximation 

JoCT-NiiCj ( nJk) (tivjjv- l^»l. 

A pair of integrals relevant to the electromagnetic power output from thin 
conducting wires are expressed by 

( 86 ) 

« 1 • Vif/ 

and 
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(87) 


the first (Oppenheimer, 1970) rests on the assumption of a longitudinal current 
that vanishes (only) at the ends of the wire and varies sinusoidally in between, 

whereas the second (Papas, 1965) is linked with a hypothetical travelling wave 

ik/Yx 

current distribution (,e ) along the wire. Both integrals involve a para- 

meter 3 = k£/2 that depends on the free space wave number k and the length i 
of the wire. Neither of the references contains any analysis bearing on the 
integrals, although their explicit evaluation is feasible. 
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The substitution x = cos 0 converts (86) to 
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and, after utilizing the differentials d(- -^) = , d(-;:r^) = 


d? 


r ,2 


(2-0 2 


to 


Initiate an integration by parts, it is eventually found that 
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■( 


in terms of the sine and cosine integrals. 

The change of variable previously employed for the integral (86) yields, 
after its application to (87), 
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whence 




Cl) 

in accordance with the particular values F) ^ (0,a) = 0 , — (O.a) = Q ; 

D dp ’ 

similarly, 
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and thence, after two integrations, 
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describes the appropriate function. Combining the expressions for and 
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Fg there obtains, finally 
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(89) 
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Newman and Magnus (1959) encounter several definite integrals in their 
study of maximal antenna gain, with the prototype 




and they offer a one-term estimate (referring to a long, though not detailed, 
elementary computation) 


I = (-1) 
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may 


be noted and its version in the case m = n is likewise available. 
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